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Abstract: We determine the structure of the hydrodynamics with conserved current, 
using the gauge/gravity duaUty of charged black- hole background. It turns out that even in 
the presence of the external electromagnetic field at the boundary, bulk Einstein equation is 
equivalent to the boundary conservation of energy momentum tensor and that of current. 
As a consequence, the thermal conductivity and electric conductivity are calculated in 
terms of the parameters of the fundamental theory. We find that Wiedermann-Franz law 
hold with Lorentz number 1/e^. 
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1. Introduction 

After the discovery of consistency of AdS /CFT [Q, ^, ^ and that of RHIC experiment on the 
viscosity /entropy-density ratio |5|, ^, much attention has been drawn to the calculational 
scheme provided by string theory. Some attempt has been made to map the entire process 
of RHIC experiment in terms of the gravity dual The way to include chemical potential 
in the theory was figured out in |9| and phases of these theories were also discussed in 
D3/D7 setup The essential feature is that the global U(l) charges in the 

boundary corresponds to the local U(l) charges in the bulk. 

The gauge/gravity duality holds for the theory with = 4 Super Yang Mill system 
which is far from the real QCD. Therefore we need to search quantities which does not 
depends on the details of the theory like rj/s. The hydrodynamics approach is especially 
useful in this respect since it is about dynamics of long wavelength and low frequency limit. 
Also, due to high temperature nature, all the fermions are decoupled. 

Recently, in a very interesting paper [^], the authors showed that the hydrodynamic 
equation in the boundary theory is equivalent to the Einstein equation in the bulk. Usually 
the dissipative parts of the current and energy momentum tensor are constructed from the 



second law by considering the entropy current [14|. The method of reference |13] has the 
merit to determine the dissipative part of the energy momentum tensor without using the 
second law. This merit was applied to other different black holes [||, |l9|, |2^. The 
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purpose of this work is to determine the dissipative part of the conserved current along 
this hne. We wih also get transport coefficients in terms of the data of the Yang-Mill 
theory. Since the dual of the particle number can be regarded as the local charge of the 
Maxwell field in the bulk, we are lead to consider the charged black holes in AdS space. 
We consider only up to the first order in the derivative expansion. It turns out that our 
expression for the current is precisely agree with that of Landau-Lifshitz [14, 15 1. The only 



difference is that the hydrodynamics transport coefficients are now determined from the 
fundamental theory. 

The rest of the paper goes as follows: in section 2, we give an outline of calculational 
procedure following |13|. In section 3, we calculate the structure of conserved current using 



the perturbative solution to the first order in derivative expansion, which is enough to 
determine both thermal and electric conductivity. In section 4, we give a summary and 
discussion. In appendix, we list the important formula. 

2. General construction 

2.1 General procedure of calculation 

In this section, we explain the general procedure of constructing solution. Our starting 
point is the Einstein Maxwell theory with a negative cosmological constant. The action is 
given by 



1 ' .5. 



levrG 



M 



+ counter terms , (2-1) 



where the first and second terms are usual Einstein-Hilbert action and cosmological con- 
stant action, and third term is the Gibbons-Hawking term, where is trace of extrinsic 
curvature Q^u- Introducing spaces with boundary, this term should be added to the action. 
In addition, the counter terms are needed to cancel divergences of the action and boundary 



stress-energy tensor. This is well studied in [16]. Taking a unit, where / = 1 and IGvrG = 1, 
the equations of motion are given by 

Rij - \gijR - Qgij - ^ (fikFj"" - \gijF^^ = (2.2) 

VjF^i = . (2.3) 

The solution we are interested in is the charged uniform black brane solution, which can 
be written as 

^ -r'^f{r)dv^ + 2dvdr + r'^{dxl + dxj + dx^) (2.4) 

fir) = 1 - — + 7^ = -6 

2^/3Q , 
F = g 2 — av A ar . 
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Here we use the Eddington-Finkelstein coordinate v = t + r^ with dr* = dr/f to make the 
metric regular at the horizon r+. It is slightly different from time coordinate 't'. In the 
boosted frame, it can be written as 

ds'^ = -r'^f{r){u^,dxi'f - 2u^dx^'dr + r'^P^^dx^'dx'' (2.5) 

F = -gt/^u^dx>' Adr , A= (^eAf - ^^^n^^ dx^ 

1 i Pi 

u = — , , u - 



\/^-Pf 

where M, Q, A'^^ and /?*'s are constants parametering solutions. It satisfies the equations 
of motion. 

Now we lift above parameters to the functions of x^. Then, (|2.5| ) is not a solution of 
the equations of motion any more. To find make it a solution, we have to add corrections 
to it. To do this systematically, first we define following tensors 

Wij = Rij + igij + ^ (^FikF^j + IgijF^^ (2.6) 
Wi = VjF''i . (2.7) 

The equation of motion is given by the vanishing of Wjj and Wj. Taking parameters 
as functions of x'^ gives nonzero Wjj and Wj. By construction these are proportional to 
derivatives of parameter functions. We will call these terms by source terms. To cancel 
this source terms, we need to add some corrections to metric and gauge fields. We consider 
the derivative expansion and do it order by order. So obtaining new solution is reduced to 
finding these correction terms. Due to the redundancy in metric and gauge field, we can 
take a gauge choice for the correction terms. Our choice for n-th correction terms are 

^^(n)2 ^ hl^dv"^ + 2K{r)dvdr + 2^-^dvdx' + [a^^ - hin(r)6i^ dx'dx^ (2.8) 
=a'i{r)dv + a'^{r)dx' . 



For metric parts, this choice is the same as [|13[. Since, for gauge field part, a"(r) does not 
contribute to field strength, the choice a"(r) = is trivial. 

Our procedure can be described as follows. First, we consider zeroth order solution. 
By construction, sfj and Sj^^ are zero, and the zeroth order correction terms are also 
zero naturally. For the first order, we take the parameters in zeroth order solution as 
functions of x^ and expand around x'^ = up to first order in derivative^. Putting these 
new metric and gauge fields into Wjj and Wi, we obtain the source terms sfj = —Wjj 
and S^P = —Wj. To make these source terms canceled, we consider correction terms in 
( |2.8| ) for n = 1, these make effects in Wjj and Wj. So Wjj = (effect from correction) — S'jj 



^We assume that 13' (0) = 



-3- 



and Wj = (effect from correction) — S^^\ from this we can find corrections. Assuming we 
know (n — 1) th order solution, we can calculate nth derivative order source terms and we 
can obtain nth order correction terms in metric and gauge fields. For general nth order, 
in the Maxwell equation, what we have to do is just solving following equations 



f(r) 



{r3a;;'(r) + AVSgQh^'ir)^' - 5(")(r) = (2. 



= 1 ir'f{r)ar{r) - ^^j;(r)| - ^("^(r) = 



For Einstein equations, the equations are 

W.. = f{r) |-i-^ar(r) - 8r'hn{r) - r\r' f {r))' (r) - (^^^'^ - S^Hr) = 
W.r = 8K{r) + {r^f{r))'h'Ar) - ^^^"'(r) + ^ (^)' " St\r) = 



rr 



hr^h'^{r))'-Si^Hr)=^ 



W,,5^^ = 2ir^Kir) + h^r) + ^^<'ir) + l{r^^ f(r)h'^{r)y - si;\r)6'^ = 

W,, - l6,,{6'^'Wki) = -^(rV(r)ar/W)' - (r) + i<5,,.(<5'='5(;)(r)) = . (2.10) 

After solving above equations nth derivative order solution can be found. 

There is an important point in this step, actually all equations in (|2.9| )and ( |2.10| ) are 
not independent. So we can find relations between some source terms 

W^+ry{r)Wr = : 5^") + r2/(r)5^") = (2.11) 

W,, + r^f{r)W„ = : S^^ + r^f{r)St'^ = . 

Thus our integration procedure is that we integrate Wr,Wi,Wvr,Wrr,Wri,Wri and Wij 
equations and consider the constraints ( 2.11| ). We anticipate fluid dynamics equations 



d^Jf^ = , dfj_T'^'^ = from the constraints. 

To get more practical guide line to our goal, now we describe way to finding correction 
terms more explicitly. Considering Wrr,Wr,Wij6^^ and traceless part of Wij, correction 



^Where ' means derivative of r coordiante. 
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terms can be obtained as follows 

hn{r) = dy^ j\''Sif{x) (2.12) 



kn{r) = I dx- 



T 

X 







Si';\x)6'^ - 2Ax^K{x) - '^Qa^J[x) - \{x^^ f{x)h'Sx))' 
■' gx x' 



al = -fdy-^^ fdx2x [st;\x) - U^.S'^S^^ 

Prom the above equations, we can find corrections scalar and tensor corrections. Using the 
result, we have to check that W^r = equation is consistent. Since the remaining equations 
VFj = and Wri = are coupled each other, it is more difficult to solve these than others. 
These equations are 

f>)'-^<.;M = S«(0. (2.13) 

M/(r)a;(r) - ^^i.(r)~j = rSf (r). (2.14) 

For the exact solutions to these equations, see appendix 
We define the chemical potential as 

// = ^o(r+)-^o(oo) . (2.15) 

One should notice that this is the chemical potential that couples with the number density, 
not the charge density. See appendix |A|. 

2.2 Stress Energy Tensor and Current 

In this subsection, we describe the prescription to calculate boundary stress-energy tensor 
and current out of bulk quantity briefly. We will follow prescription in for stress energy 
tensor and calculate the boundary current. In order to calculate these tensors, we need to 
decompose our metric. Our metric decomposition is a well-known ADM decomposition 

ds"^ = -i^u{dx^' + V^'dr){dx'' + V''dr) + N^dr^ , (2.16) 

where is defined by ^^"Vy, and the 7^^^ is the inverse of ^^y. Under this decomposition, 
we can obtain the stress energy tensor from the action ( |2.1D by variation of boundary metric 
^fj^y. The stress energy tensor is given by 



2 6 S I 1 
.f,^ = ..... r^—=—^ = lim r'^[-2{@^,y - 67^,, + 87^,, + -G^„y) 



1 — >oo 



where the last two terms came from counter lagrangian which has been studied in |16| 
and subscript "c?" means that we have to impose equation of motion. In our metric 
decomposition, the expression of extrinsic curvature Q^y is 

Q^.u = ^[l^.-D^V,-D,V^] . (2.18) 
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Another important physical quantity is the boundary current. The current can be obtained 
by same method of the stress energy tensor. The current is 



= Hm r^^^ = Um r^^F^^ 



(2.19) 

-7 5A^ r-^cc 



= lim r^^[j>^\A\ - dxAr) - V^^^^dxA^ - d^Ax)] 
where An is the gauge field which is projected to the boundary. 



3. Perturbative solution and fluid dynamics 
3.1 Perturbative solution 

As we discussed in previous section, to obtain higher order derivative solution, we should 
take the parameters of the solution as functions of boundary coordinates. Let us rewrite 
the zeroth order solution with functions. The zeroth order solution is 

gfjdx^dx-^ = -r'^f{r)u^u^dx^'dx'' - 2u^dx^'dr + r'^P^^dx^'dx" (3.1) 
= [eA^^\x) - ^^^^u,{x)^ dx^ . 

This solution satisfies the Einstein equation and the Maxwell equation in zeroth order. 
Notice that we need to distinguish the R-charge coupling g and the external electric charge 
coupling e. Using the gauge field part and the definition of chemical potential ( 2.15| ), we 
can read chemical potential as 

M:,0 = 4^ ^ (3,2) 

Putting this zeroth order solution into the general formula of boundary tensors, (|2.17| ) 
and ( gloD , we can get the zeroth order boundary (particle number) current and energy 
momentum tensor, 

^ .= nu'' (3.3) 

lUj g 

= 2M(ry'^^ + 4n^ii'^) . 

These are well-known perfect fluid current and energy momentum tensor. 

Next step is to obtain first order corrections. To construct first order correction, we 
regard parameter M,Q,A^^^ and /3j as functions of x^ in ( |2.5[ ) and expand these parameter 
functions around x^ = 0, then this expansion gives source terms sfj{r,x) and sf\r,x). 
After a little working, we can find first order source terms. The source terms S^P and S'jj 
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are given as follows 



^vv (^) 


z 




_ di(5i 

T 




~ 




= 




3 







3 ^ 5 



(3.4) 



2r / j,o g 



3 
2 



4'^ (0 = 5^ {d,Q + QdiPi) (3.5) 
5«(r) = 

si'\r)=9(-^id^Q + Qd,Pi)---F:A , 
\ r-^ r g I 

where F!^f^ = djjAf^^ — diAf^^ is the external field strength tensor. Using above result, one 
can easily calculate parts of correction terms, 

/i(i)(r) = , = , fc(i)(r) = lr^dil3' (3.6) 

(n) ^ „ ^ . 2 



where definition of a{r) and it's asymptotic expression are given by 



-M = 3y^*p^y^^<i..^«<;-j^ + ^-^-^ + o^^-j . (3.7) 

There is another constraint W^r = which must be checked by the above information, and 
it is very easy to check that the above correction terms satisfy the constraint. 

Though we can write explicit expressions of ji{r) and ai(r) from appendix it is 
too long to write in the main text and we only need the asymptotic behaviors of those 
solutions. In order to write covariant form, we decompose first order vector terms as 

a,(r) = ap{r)d,Pi + aQ{r){diQ + Qd,(5i) + (r)F;f , (3.8) 
ji{r) = jp{r)^,P^ + jqirMQ + Q9,A) + jF(r)F;f. (3.9) 

With this decomposition, the gauge field can be written as 

AW(r) = ap{r)u'^d,u, + aQ{r)u- F^f + ap{r)u^Ftf, (3.10) 
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where F^'' = d\ (Qu^) — dy {Qu\). The form of metric is 



ds^ = —r^ f {r)u^Uydx^ dx'^ — 2u^dx^dr + r"^ P^^dx^ dx'^ + 



2r 

—u^Uydxu^ (3.11) 



dx^'dx"" 



where we took the covariant expression and the definitions of 0"^,^ is 

a^' ^\p^'''P'^{daUp + dpUo,)-^P^"'dau'' . (3.12) 
The chemical potential in 1st order calculation is 

■ (3.13) 

where 5r^^ is the first order correction of the horizon. As a final step for completing our 
consideration, we have to find (5r^\ so we need to examine the first order correction of 
horizon. In general, the horizon is defined by following equation 

giJdj{r-rH)dj{r-rH) = ^ ■ (3.14) 

In zeroth order, we already know the horizon r// = r+. As a check, zeroth order equation 
is just g'^'^ = 0, then it's solution is the right value r+. In the first order, the equation is 



g^' - Ig'^d^m = 0, (3.15) 

ion to the 

keeps same expression, 



but there is no correction to the solution and Sr^^ vanishes. So the chemical potential 



VSgQis 
r'^{x) 



where we have to notice that Q and r^. are not constants. 

3.2 Boundary tensors and fluid dynamics 

Putting our source terms into the equations ( |2.11 ), they give constraints 



+ ^MdiPi = (3.17) 

diM + 4Ma„A = VSQ-F^. 

9 

d^Q + QdiJ3i = . 



■pcxt 
vi 
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Figure 1: The coefficient of tliermal conductivity or The electrical conductivity 



we can rewrite these in covariant form, 



(3.18) 



Since there is a derivative, these are exact first order (non)conservation equations. The right 
hand side is energy and momentum inflow sourced by the external fields. When external 
fields are absent, we obtain the usual fluid dynamics system with conserved current. We 
expect that the relation 



(3.19) 



holds to all order. 

We know the asymptotic forms of metric and gauge fields. This information is enough 
to calculate the first order boundary stress-energy tensor and the current. Using our first 



order metric and equation (2.17), we can find the first order stress-energy tensor 

T^j. = 2M{r]^^ + 4:Uf,Uu) - 2r\a^u := P{r}^u + Au^Uu) - 2T](j^y , (3.20) 



from which we can read off the viscosity r] = r^. By restoring l/(167rG) = N^/Stt'^, which 
was set to be 1, we get 



''"87r2^+" 8 I2 + V4 eV^^rT, 



(3.21) 
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where we have used r+ = ^ ^1 + y 1 + -^{fi/ girT)'^ j. In addition, we can also calculate 
the first order boundary current from ( 2.1S| ). The current is 

\ r\ r. 

(^-2^Q^^-^ + ^r+^n^F-V|, (3.23) 

where the j^(r+), jQ(r+) and jpir-^) are values of each function at the horizon, which can 
be read from ( B.16 ) as follows 

(r+) _ 2 {2rl + Q^) jg (r+) _ Q jp (r+) _ e V^Q 





- -^(0) 


•^(1) 












(3.24) 



To find physical quantities, we need to change the above equation into familiar form. Using 
(3.18) and the temperature of the black hole, we can obtain the final expression of the 
current which contains dissipative current and contribution from the external field. After 
some algebra as summarized in the appendix the expression turns out to be simplified 
and given by 

^S. = -l?Mt^""^'f + l^"^=^"V . (3.25) 

:= -KP'^'^d^^ + -u^F^"*a'^ , (3.26) 
T e 

where 1/e in cr/e is inserted to read off the electric conductivity from the electric current 
eJ rather than the number current J. See appendix 

Now the coefficient of thermal conductivity k and the electrical conductivity are given 

by 

""-^^^ VMS • (3-27) 

To see that the second one is actually the electric conductivity, we notice that = 
au^F^^^ becomes in non-relativistic limit Jcxt = <y{E + vxB) with the a given above. It 
is easy to show that k can be also written as 

- - ^+ (2 - a)2 ^.^^ ^ ^ ^2/^6 (3 28) 



4g2 (l + a)2^ 



so that our result agrees precisely with that of |22], which was obtained from the Kubo 
formula. Our result of conductivity goes to the known result in |24, 25] for zero charge case 
^. The relation between the thermal and heat conductivity, so called Wiedermann- Franz 
law: 

K = aT/e^, (3.29) 



^The result is difFerent by a constant factor of tt. 
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is manifest in our case with Lorentz number l/e^. Considering the comphcation of the 
definition of temperature in the charged black hole and also the difference in the origin 
of two terms, it is rather surprising to have such similar and simple recombination such 
that Ta and k is proportional to each other. Using 2M = r^(l + fi"^ /{3g'^r'^)) and r+ = 
7rT(l + Y^l + 2{fi/g7rT)'^ /3)/2, we can express k/T^ and a/T as a function of fi/gT only: 



rp2 g2 



i-3v/i + i(^)^)' 



Fig. H shows this as a function of fJ./gT. 

We can also calculate the thermal conductivity, 



Notice that = and is 1 and N(./2Nf for R-charge and Baryon charge respec- 
tively pi. 



4. Discussion 



We use gauge/gravity duality to determine the structure of the fluid dynamics in the pres- 
ence of the conserved current and calculated thermal conductivity as well as the electrical 
conductivity in the presence of the the external electric field. Since the dual of the par- 
ticle number can regarded as the local charge of the Maxwell field in the bulk, we used 
the charged black holes in AdS space. For our purpose, the first order in the derivative 
expansion. 

While the determination of the current especially the dissipative part is of highly 
interesting, going second order or higher order is less interesting from the experimental 
point of view: what RHIC experiment discovered is that we should neglect the dissipation 
part almost completely. This perfect liquid behavior is a hall mark of the RHIC experiment. 
However, LHC experiment may show different behavior due to its much higher collision 
energy scale. Therefore calculating the transport coefficients may be of some importance 
for future experiment. It would be also interesting to calculate the similar quantities for 
other ads space for the application to the solid state physics or M2, M5 brane theories. It 
would be also interesting to see if this method can help to settle problems of fluid mechanics. 

m. 
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Note added 

At the end stage of the this work, we received two papers ^] whose contents are partly 
overlapping with present paper. Therein, external fields are not included but T^u and 
are calculated up to second order. 



A. Number current v.s charge current: Fixing g factors. 

For the lagrangian 



1 ^2 



the equation of motion is 



L = -j^F^, + J,A'^ (A.l) 



^d^'F^, = J, (A.2) 



g 

If we rewrite the action in terms of the rescaled variable A^/g := A'^, 

L = -^F'l,+gJ^A'^', d^Fl,=gJ,:=Jl. (A.3) 

therefore if j'' = n is the number density, Jq = gn is the charge density and is the 
number current. Also the electromagnetic potential is A'^ = A^/g. 

The number current density = 2\/3Qui^/g so number density is related to the charge 
parameter Q by 



n 



2V3Q/g. (A.4) 



The chemical potential for the electric charge density is /u' = A'{qo) — A'{0) = ^/SQ/r'^. 
From the A.J = A'. J' or fin = h'Jq and Jq = gn, we get /i = gfi' . Then the chemical 
potential fi that couples to the number density is /i = Aq{oo) — ^o(O) as was stated in the 
main text ( p. 151 ). Similar relation hold for electric charge as well as the R-charge. 

B. The exact expression for of ji{r) and ai(r) 

In this appendix we will present exact expressions of ji{r) and aj(r) which can be obtained 
from coupled differential equations ( |2.13D and ( |2.14| ). Integrating ( |2.14| ) from r = r+ to 
r = cxD, we get 

r^fir)a',ir) - iVSgQ - ^-^'^ = £ dxxSi{x), (B.l) 



where we used the fact that /(r) is zero and ai{r) is not singular at r = r_|_. Using ( 2.12| ) 
to eliminate a^(r) in ( |B.1| ), we get a second order differential equation of ji{r) 

- -/^r) - ^f(r) = Ur), (B.2) 

(B.3) 
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where 

= -3r8.ft - (2v/3Q=!iij:i - V3 ( - - —] (d.Q + QOM 

+ (r - r+) ) . (B.4) 

Two linearly independent homogeneous solutions of this equation are 

JiTiW = rV(r) (B.5) 

poo ^3^„ poo J 

A particular solution to (|B.2| ) is found by using method of variation of parameters 

jp{r) = bi(r)jH^{r) + b2{r)jH2{r), (B.7) 

where 

b,ir) = - r d.i^^^^^, (B.8) 

b,{r) = r^d,(3, + dx {iElk^^l^ + Sx^a.ft) . (B.9) 

Note that Sx'^d^Pi term cancels the divergence of the integral. Now we can write the general 
solution to 



Ur) = jp{r) + CijuAr) + C23H,{r)- (B.IO) 
To determine Ci and C2 we need to find the asymptotic behavior of ii{r) 

J?ir) = C,r^ + r^S^A + + Q (1) . (B.ll) 

The coefficient Ci multiplies a non-normalizable metric deformation, and so is forced to 
zero by our choice of boundary conditions. The other integration constant C2 leads to a 

(n) 

nonzero value for Tq^. We can remove this ambiguity by demanding that u^Tj^u = 0, thus 
C2 can be set to zero. In summary, 

dy 



3i{r) = -r f{r) / dxxf{x)Ci{x) 



+rV(r) ^57^) (^'^-/^^ + /"^^ [xfi^)Ux) + Sx^d.Pi]^ .(B.12) 



Having the expression of ji{r), ai{r) is obtained by integrating ( |2.13| ) 



ai{r) = -g dr 



00 ^3 



\/3Q 



X" X ( i[{x) 



2 V x3 



Sri{x) 



(B.13) 
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where we chose the gauge to make aj(r) vanishes at infinity. 

Since we are interested in boundary physics, we present asymptotic behaviors of these 
functions here: 



ji{r) « r3a,A-^^-F-t + 



2V3Q e 

r VI 

5r g 



2 \ 2r_|_ 



r+Qe ext 



(B.14) 



a,(r) « -eF:f + ^ I VSgQ^-^^ + ^ (9,Q + Q9,ft) - ^eF^f* ) . (B.15) 



2r+ 



Still we have undetermined constant, jj(r+). Taking r — > r+ limit to ( |B.12|) yields an 
equation of jj(r+). After some calculations we get 



(r+) _ (2r4 + r^rl + rt) d,f3i - Q (diQ + Qd^Pi) - ^/Sr^ g|F; 



ext 
vi 



6r+Q2 _ 2rl (-2r^ + r?. + ri) 



8Mr3 

From the definition 

ji{r) = jf3{r)^,P^+3Q{r){^^Q + + ji.(r)F:f , 

we can read off the result ( p.24[ ) . 



(B.16) 



(B.17) 



C. Current expression and useful relations 



To obtain the current expression of ( 3.25 ), we massage ( 3.23 ) and get following form 



(1) 



2 

4Mr+ 



(C.l) 



We need to change first term into derivatives of physical quantities. For this, we use more 



useful form of ( 3.18| ), 



AM g 



AM 



-P^"'^yM . 



(C.2) 



Using this equation and other relations between physical quantities, one can get our result 
(3.25). We write down the useful relations for this algebra as 



Q 



V3g' 



2 2,4 A* ^+ 

r_|_r_ + r_ — 



[2r'\ — r^r^ — r 



Qrl - fiyg^ 3rl - 2M 



67rr_i 



27rri 



(C.3) 
(C.4) 
(C.5) 
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